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A  ROBUST  ADAPTIVE  ARRAY  STRUCTURE 
USING  THE  SOFT  CONSTRAINED  LMS  ALGORITHM 


Gary  F.  Hatke  and  Stuart  C.  Schwartz 

Department  of  Electrical  Engineering 
Princeton  University,  Princeton  N.J.  08540 


ABSTRACT 

This  paper  presents  a  new  algorithm  for  adaptive  beamforming  which  is  robust  to  both 
desired  signal  misalignment  and  array  geometry  imperfections.  It  extends  the  work  done  by  Jab- 
Ion  and  Widrow  [ij  by  utilizing  a  priori  knowledge  of  the  desired  signal  direction  to  generate 
new  soft  constraints  on  the  pattern  of  the  array. 


1.  INTRODUCTION 

Adaptive  beamformers  can  be  used,  with  minimal  a  priori  information,  to  greatly  improve 
the  signal- to- noise  ratio  of  a  signal  in  the  prescence  of  hostile  interference  (jammers).  However, 
when  the  actual  signal  angle  of  arrival  does  not  coincide  with  the  estimated  angle  of  arrival,  or 
when  the  actual  geometry  of  the  array  differs  from  the  geometry  assumed  by  the  array  designer, 
then  the  performance  of  the  array  may  be  degraded  due  to  adaptation.  This  phenomena  is  due 
to  the  fact  that  in  the  absence  of  a  pilot  signal,  the  usual  adaptive  array  attempts  to  minimize 
total  output  power,  regardless  of  the  type  of  signal  present.  While  this  is  desirable  if  the 
received  signal  is  mostly  jammer  and  receiver  noise,  when  the  desired  signal  component  becomes 
large,  the  receiver  attempts  to  reject  it,  thus  significantly  degrading  performance  [2,3]. 

The  most  common  way  of  avoiding  the  problem  of  rejecting  the  signal  along  with  the  noise 
is  to  impose  some  sort  of  linear  constraints  on  the  ways  in  which  the  array  is  allowed  to  adapt. 
Many  methods  have  been  proposed  to  implement  these  types  of  algorithms,  but  the  most  com¬ 
mon  are  the  Frost  algorithm  and  the  Generalized  Sidelobe  Canceller  (GSC)  structure  (4,5).  Both 
of  these  types  of  arrays  impose  constraints  on  the  array  patterns  which  they  may  generate,  usu¬ 
ally  either  fixed  direction  gain  constraints,  or  constraints  on  the  derivative  of  gain  vs.  angle  at  a 
fixed  look  direction. 

The  array  processor  proposed  in  this  study  is  designed  to  give  the  maximum  number  of 
degrees  of  freedom  for  the  processor  to  null  out  interferers,  while  being  very  robust  with  respect 
to  both  array  geometry  and  signal  direction.  This  is  achieved  by  use  of  the  available  a  priori 
knowledge  of  the  signal  direction  to  generate  soft  constraints  on  the  processor  weights. 
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2.  SYSTEM  MODEL. 

For  the  remainder  of  the  paper,  the  structure  of  the  adaptive  array  under  consideration 
will  be  that  of  a  zeroth  order  GSC.  This  structure  has  been  shown  to  converge  to  the  same 
steady  state  solution  as  the  Frost  constrained  adaptation  algorithm  with  one  directional  con¬ 
straint.  We  will  be  interested  in  the  case  of  narrowband  systems,  with  sources  emitting  only 
plane  wave  signals  (i.e.  far  field  approximation).  We  will  assume  one  jammer  with  unknown 
direction,  and  a  signal  whose  direction  is  random,  but  not  uniformly  distributed  over  the  entire 
(— rr.Tr)  angle.  We  also  assume  that  the  antenna  array  is  linear,  with  uniform  X/2  spacing 
between  antennas.  However,  the  gam  of  each  antenna  is  assumed  to  be  a  complex  gaussian  ran¬ 
dom  variable  with  inphase  and  quadrature  parts  both  distributed  as  N(l,(7^j),  all  uncorrelated 
with  each  other.  This  model  is  also  appropriate  for  an  array  where  the  antenna  placements  are 
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random,  but  close  to  linear. 

It  should  be  noted  that  the  assumptions  of  a  uniform  linear  array  and  one  jammer  are 
made  purely  for  mathematical  expediency,  and  are  not  necessary  for  the  results  of  this  paper  to 
hold. 

The  sidelobc  cancellation  matrix  B  is  designed  to  be  an  N  x  N-1  matrix  of  rank  N-1,  such 

that 

Bs‘='0  (1) 

where's^  is  the  vector  of  delays  seen  for  a  plane  wave  signal  coming  from  the  constraint  direc¬ 
tion.  Normally,  the  array  is  assumed  to  be  pre —steered  so  that  the  constraint  direction  is  broad¬ 
side  to  the  array,  and  thus  : 

(2) 

The  output  of  the  GSC,  using  our  assumptions  on  the  gain  of  the  antenna  elements,  is 
given  by  the  expression; 

T'GT-w'BG'x  (3) 

where 'jr  is  the  received  vector,  and  w  is  the  N-1  dimensional  vector  of  weights  used  in  weight¬ 
ing  the  outputs  of  the  B  matrix. 

G  is  the  diagonal  matrix  of  the  antenna  gains,  and  can  be  written  as: 

and  Api  is  a  N(0,<7^j)  complex  random  vector. 

We  now  define  the  output  of  the  beamforming  matrix  as  TT,  and  look  at  the  optimization 
problem  of  minimizing  the  total  output  power  of  the  GSC,  which  translates  to  the  problem; 

£:[  iT'GT-^'ir  p]  (4) 

The  solution  to  this  problem  is  the  well  known  Weiner  Filter  solution.  The  optimum 
weight  vector  w  is  given  by  the  equation: 

w  (5) 


If  we  look  at  the  case  of  a  single  jammer,  a  desired  signal,  and  uncorrelated  receiver  noise, 
then  the  matrix  R„b  looks  like: 

R«»  =B  [<72J-(-GR„G^ -f-^T^^GTTG^B^  (6) 

Where  the  vectorsV  and  j  are  vectors  of  delays  seen  at  the  antennas  for  plane  wave  sig¬ 
nals  coming  from  the  desired  signal  and  jamming  directions  respectively.  The  matrix  R„  is 
defined  as: 


R„=o^.E[T^\ 

And  the  vector  is  of  the  form; 

r„i=BGR„^GtT 


(7) 


(8) 


or 


Where  R„  =  E  [  xlT 


□ 


And  andlc"*  denote  complex  conjugate  transpose  and  complex  conjugate  of  IT  respec¬ 


tively. 


3.  ROBUST  ADAPTATION  ALGORITHM 
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If  the  signal  direction  vector  s  is  not  exactly  equal  to  1,  or  if  the  variance  of  the  antenna  — 

gains  is  not  very  small,  the  array  pattern  generated  by  (4)  will  attempt  to  null  out  a  strong 
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desired  signal.  To  deal  with  the  problem  of  random  antenna  gains,  Jablon  and  Widrow  Ilj  have 
proposed  using  a  soft  constraint  type  of  optimization.  A  suitable  penalty  function  is  chosen,  and 
the  ojDtimization  problem  is  then  changed  to  minimizing  the  original  functional,  (  in  this  case 
E\  1  1‘G3c’  —  w'TT  1*  ]  )  plus  some  penalty  constant  times  the  penalty  function.  If  the  penalty 
constant  is  chosen  to  be  large  enough,  then  the  optimization  problem  will  have  a  solution  very 
close  to  that  of  a  hard  constraint  problem,  where  the  constraint  is  defined  as  the  penalty  func¬ 
tion  being  forced  to  take  on  the  value  zero. 

The  penalty  function  chosen  in  [1)  attempted  to  minimize  the  leakage  of  the  signal  com¬ 
ponent  into  the  cancellation  branch  due  to  random  antenna  gains.  Here,  we  extend  this  tech¬ 
nique  to  the  case  of  both  unknown  antenna  gains  and  unknown  signal  directions,  and  derive  a 
new  LMS  algorithm  for  this  case.  VVe  then  show  how  the  results  can  be  extended  to  give  a  rea¬ 
sonably  efficient  algorithm  by  choosing  the  correct  type  of  beamforming  matrix. 

We  will  define  our  penalty  function  p(w)  as  one  which  is  linear  in  the  amount  of  power 
that  is  transmitted  through  the  cancellation  branch  of  the  GSC  due  to  the  signal  component. 
The  leakage  of  the  signal  vector  into  the  cancellation  branch  is  given  by  the  expression; 

L(s*)  =  w'BGT  (9) 

We  will  take  the  expected  value  of  the  power  of  the  leakage  signal  as  our  penalty  function. 
Thus,  p(w)  becomes; 

p(w)  =£;[  w‘B(  GiT^G^  +  1  (10) 

Bringing  the  ensemble  expectation  inside  and  utilizing  the  assumptions  that  the  errors  in 
the  antenna  gains  are  independent  of  each  other,  and  of  the  possible  signal  directions,  and  also 
that  the  variance  of  the  components  of  the  signal  vector  are  all  equal  to  some  constant,  c^,,  ,  we 
get  the  simplified  result; 

p(^)=^‘B(<z2,o2„H-R„)Bt^*  (11) 

And  now  the  soft  constrained  optimization  problem  becomes; 

T‘R„T-2Re[w‘BR„Tl-l-w‘B(R„-l-c(c72,(72^,I-l-R„)]B^w*  (12) 

At  this  point,  we  note  that  the  term  is  quadratic  in  the  vector  w ,  and  thus  any  minima  of 
this  expression  will  be  global.  The  matrix  R„  is  the  spatial  autocorrelation  matrix  for  the 
desired  signal.  This  matrix  is  assumed  to  be  known  o  priori,  but  as  we  will  show  later,  any  rea¬ 
sonable  estimate  will  suffice  for  the  proposed  algorithm. 

We  now  take  the  instantaneous  approximation  to  the  gradient  vector  of  our  error  surface 
using  the  usual  LMS  methods  and  subtract  some  multiple  of  it  from  our  present  solution  for  w. 
This  method  gives  us  a  new  adaptive  algorithm  for  the  tap  weights; 

=  (  I -2^( -b^B[<72„<72^I-|-R„)BMwt+2^'TBT*  (13) 

Or,  after  re-arranging  terms; 

=  [  I  -  -  iiBR„B^  ]  w ^  -f  2pyj^\  (14) 

Where  the  term  y*  is  the  output  of  the  array  at  the  kth  iteration.  This  algorithm  can  be 
shown  to  converge  for  positive  values  of  p  such  that: 

P  <  - i -  (15) 

fr!R„„+-^B(o2^a2„I+R„)BM 

It  should  be  noted  that  the  matrix  BR^B^  is  Toeplitz,  non-negative  definite,  and  hermi- 
tian.  Thus,  the  matrix  is  diagonalizable  by  a  unitary  matrix.  Let  this  matrix  be  denoted  by  P. 
Then,  we  can  write  (14)  as; 
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wt+i  =  [  I  -;;PA„P‘  jw*  +  2idyk'u\  (16) 

Where  the  matrix  A,,  is  the  diagonal  matrix  of  eigenvalues  of  the  matrix  BRj^B^._VVe  now 
change  coordinates  by  using  the  operator  P,  and  define  a  new  beamforming  matrix  B  by  the 
expression: 

B=PB  (17) 

The  matrix  B  has  the  required  property  of  a  beamforming  matrix,  i.e.  equation  (1)  is  true  for  B 
if  it  is  true  for  the  matrix  B.  Using  this  new  beamforming  matrix,  the  update  equations  now 
become: 

Wi+i  =  [  B^  ]w*  %  (18) 

Where  the  vectors  w  andu*' are  defined  in  the  new  coordinate  system. 

4.  RESULTS 

The  new  algorithm  requires  the  knowledge  of  the  signal  spatial  autocorrelation  matrix  R,, 
a  priori.  In  practice,  this  can  be  approximated  by  assuming  that  the  signal  will  definitely  fall 
within  the  main  beam  of  the  array  (or  else  the  array  would  be  mis-steered),  and  then  taking  the 
maximum  variance  distribution  within  this  arc.  Using  this  method,  the  spatial  distribution  of 
the  desired  signal  is  assumed  constant  within  some  arc  symmetric  about  broadside,  and  zero 
elsewhere.  Since  this  is  an  even  symmetric  function,  the  matrix  R„  has  all  real  entries. 

The  quiescent  pattern  generated  by  the  new  adaptation  algorithm  was  calculated  analyti¬ 
cally  for  the  case  of  an  11  element  linear  array  in  the  presence  of  one  signal  and  one  jammer 
The  SINR  was  calculated  as  a  function  of  jammer  angle  of  arrival  for  a  fixed  signal  location, 
and  then  again  as  a  function  of  signal  angle  of  arrival  with  a  fixed  jammer  location.  The  results 
show  that  the  loss  in  SINR  due  to  lack  of  exact  knowledge  of  the  desired  signal  direction  is  very 
small  as  long  as  the  jammer  lies  outside  the  main  lobe  of  the  array  pattern,  and  the  actual  sig¬ 
nal  vector  lies  inside  the  main  lobe.  This  is  due  to  the  fact  that  the  soft  constraint  "freezes"  the 
pattern  of  the  array  in  the  main  lobe  region,  while  allowing  it  to  adapt  in  other  regions. 

A  comparison  of  the  output  SINR  of  the  new  soft  constrained  algorithm  and  a  first  order 
constrained  GSC  shows  that  the  new  algorithm  is  much  less  sensitive  to  errors  in  desired  signal 
direction  of  arrival.  In  addition,  trade-offs  between  ability  to  null  jammers  close  to  broadside 
and  robustness  to  signal  misalignment  can  be  made  through  changes  in  R„. 

5.  CONCLUSION 

This  paper  develops  a  new  algorithm  for  robust  beamforming  in  the  prescence  of  uncertain 
signal  directions  and  array  geometries.  The  new  algorithm  incorporates  available  a  priori  infor¬ 
mation  ii.to  a  penalty  function  which  imposes  soft  constraints  on  the  adaptation  of  the  array. 
The  proposed  algorithm  is  shown  to  be  better  than  linearly  constrained  adaptive  arrays  at  han¬ 
dling  off  center  desired  signals. 
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Fig.  1:  Output  SNR  vs.  Jammer  Angle 


Results  for  a  linear  11  element  array 
Width  of  signal  spatial  autocorrelation 
matrix  for  robust  algorithm  is  in  radians. 

Fig.  1-3:  Output  SNR  for  adaptive  processors 
using  infinitely  slow  adaptation. 

Fig.  4-5:  Array  patterns  due  to  simulated  gaussian 
noise,  jammer,  and  signal  environment. 


Fig.  2;  Output  SNR  vs.  Jammer  Angle.  No  Misalignment 
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Fig.  3:  Output  SNR  vs.  Signal  Angle,  Jammer  at  .27  Rads. 


Fig.  4:  Simulated  pattern  due  to  Robust  algorithm 


Fig.  S:  Simulated  pattern  due  to  0th  order  CSC. 
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